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Motionless domains walls representing heteroclinic temporal or spatial orbits typically exist only 
for very specific parameters. This report introduces a novel mechanism for stabilizing temporal 
domain walls away from the Maxwell point opening up new possibilities to encode information in 
dynamical systems. It is based on anti-periodic regimes in a delayed system close to a bistable 
situation, leading to a cancellation of the average drift velocity. The results are demonstrated in a 
normal form model and experimentally in a laser with optical injection and delayed feedback. 


Spatially extended nonlinear systems often admit mul¬ 
tiple coexisting stable states, and the dynamical proper¬ 
ties of the fronts connecting them are fundamental in the 
understanding of pattern formation and localized struc¬ 
tures (LS). These dissipative objects [1-11] occur in many 
natural and laboratory systems and are characterized by 
a correlation range much shorter than the size of the sys¬ 
tem, making them individually addressable. Localized 
structures due to a front between a homogeneous and a 
cellular state are intrinsically stable [12-14], as the pat¬ 
tern oscillations stabilize the position of the front against 
small perturbations. In contrast, fronts between homo¬ 
geneous states are typically stable only for very specific 
parameters and lead otherwise to a coarsening. 

In this letter, we discuss a robust mechanism for the 
stabilization of fronts between homogeneous states based 
on a generic property of time delayed systems having no 
equivalent in spatially extended systems. We show that 
structures occurring within a delay time can be stabi¬ 
lized by a temporal inversion after each time delay, thus 
allowing the stabilization of self-localized temporal do¬ 
mains in temporal slices modulo twice the delay. The 
case of equispaced domain walls reduces to the so-called 
square-waves [15, 16]. This has not only relevance for the 
fundamentals of the dynamics of delayed system, but also 
meets recent interest in temporal LS realized in photon¬ 
ics [17-23]. Due to their intrinsically fast dynamics, the 
possibility to use LS as bits for information processing 
was addressed early in nonlinear optics [24-26] . Interest¬ 
ing results were achieved for spatial solitons in semicon¬ 
ductor microcavities [27-29], although it turned out that 
spatial disorder limits parallelism and control [30, 31], 
motivating further studies into temporal LS. 

As indicated, in the simple case of a unidimensional 
bistable system with a single dynamical variable 'ip (x, £), 
the stable coexistence between two homogeneous phases 
is merely achieved for a single parameters’ value, the so- 
called Maxwell point. Here a domain wall separating the 
two states would be motionless. Yet, such a regime pos¬ 
sesses little experimental significance since any deviation 
of the control parameter or any symmetry breaking effect 
implies that one of the two bistable phases will eventu¬ 


ally invade the other in a way reminiscent of nucleation 
bubbles in first order phase transitions. 

In recent years, building on the strong analogies be¬ 
tween spatial and delayed systems [32-34] , a similar sym¬ 
metry breaking induced coarsening dynamics was shown 
to occur in delayed bistable systems [35, 36]. The ability 
to control the motion of these walls would have signifi¬ 
cant implications as for instance to encode and process 
information with a fast nonlinear delayed system. Mo¬ 
tivated by this idea, the pinning of domain walls was 
recently demonstrated via an external temporal modula¬ 
tion [37]. Departing from the analogies with spatial sys¬ 
tems, we demonstrate in contrast a stabilization mecha¬ 
nism based not upon a fast active modulation, but upon 
a slow, self-induced dynamics. We envision the use of a 
lesser known property of delayed systems: their ability to 
generate anti-periodic output, i.e. temporal traces get¬ 
ting inverted after each time delay r, thereby inducing 
an effective periodicity 2 r. 

While this effect was studied theoretically in [15, 16, 
38] and square-waves were demonstrated experimentally 
in several optical [39, 40] and opto-electronic systems 
[41], we demonstrate in this letter how this generic prop¬ 
erty of anti-periodicity can be harnessed to create robust 
motionless domain walls and prevent the coarsening phe¬ 
nomenon. This allows to store information even far from 
the Maxwell point and/or in the absence of bistability. 
Such an idea has no equivalent in bona fide spatial sys¬ 
tems since the anti-periodicity would actually correspond 
to a space defined over a Mobius strip. We evidence ex¬ 
perimentally and theoretically stable domains in an in¬ 
jected semiconductor laser with delayed feedback, that 
are insensitive to symmetry breaking and exist beyond 
the bistability. 

We base our analysis on the normal form of the imper¬ 
fect pitchfork bifurcation modified by the inclusion of a 
linear delayed contribution 

^ + /3$ 2 - V 3 + ri'i> (t (1) 

In Eq. (1), the parameter fi controls the bistability, /3 
represents the symmetry breaking and 77 is the ampli- 
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tude of the time delayed feedback. The temporal scale 
is normalized by the time delay r. We study the limit 
of long delays r —>• oc, and define a smallness param¬ 
eter e = 1/t making apparent the singular nature of 
Eq. ( 1 ). Delay Differential Equations (DDEs) possess 
in the long delay limit an eigenvalue spectrum that can 
be divided in two parts, see for instance [42] and ref¬ 
erences therein. A quasi-continuous branch stems from 
the influence of the delayed contribution while a discrete 
spectrum is generated by the instantaneous linear terms. 
We are interested in the regimes where a portion of the 
quasi-continuous branch may become unstable giving rise 
to smooth dynamics. We set p < 0 ensuring the stability 
of the discrete spectrum. When £ —)► 0, the left hand side 
of Eq. (1) can be assumed, as a first approximation, to 
vanish. Such an approach is enlightening as one finds a 
functional mapping governing the evolution of the small 
deviations from the trivial solution 4/ = 0 as 


T {t) = -ftT (t - 1 ) ( 2 ) 

with k = Due to the infinite dimensionality of 

DDEs, an initial condition must be given as a function 
defined over an interval equal to the delay. If |ft| < 1, any 
initial condition slowly decays from one round-trip to the 
next and the steady state T (£) = 0 is asymptotically sta¬ 
ble. Yet, for |ft| > 1, the trivial solution may bifurcate 
via two widely different scenarios. If ft < —1, the tem¬ 
poral profile evolves regularly from one round-trip to the 
next while for ft > 1 , the profile gets inverted at each 
round-trip signaling the onset of a period-2 (P2) regime. 

The slowly evolving dynamics of the temporal profile 
can be better understood and visualized via the spatio- 
temporal equivalence between delayed and spatially ex¬ 
tended systems [16, 32-34, 42]. We define the deviation 
from the two bifurcation points p = kt] with ft = =b 1 , 
as p — a + £ 2 ai and r] = aft + b\£ 2 and, as detailed in 
[42], inserting a multiple time scales expansion in Eq. (1) 
for both the temporal derivative and the delayed term, 
and defining T = sip + O (£ 2 ), one obtains the following 
equivalent partial differential equation (PDE) as a third 
order solvability condition 


dip 

dn 


, , ^ , 2 , V’ 3 , 

pip H —ip H-b 

a a 


1 d 2 ip 
2 a 2 dx 2 ’ 


( 3 ) 


with p = — (ai — ftfri) /a. Although strictly valid in the 
vicinity of the bifurcation points, such normal forms are 
known to have a wider domain of validity, see the dis¬ 
cussion before Eq. ( 8 a) in [16]. In Eq. (3), the spatial 
coordinate ( x ) must be understood as a local time co¬ 
ordinate within the round-trip while the slow time ( n ) 
represents the evolution of the temporal profile from one 
round-trip to the next. Besides, we factored out a drift 
velocity defined as v = a -1 as in [23, 33], representing the 
small deviation of the period with respect the time de¬ 
lay. Equation (3) must be complemented by a boundary 
condition that reads 


( 4 ) 



Figure 1. (color online). Illustration of dynamics in space-time 
diagrams. These space-time diagrams are folded over a time T = 
t T |?7 1 1 • a) Temporal pattern imposed as an initial condition. 
The PI regime where k = — 1 leads to the coarsening dynamics 
of the temporal time trace depicted in b) and c) towards the high 
and low states when /3 = — 10 -2 b) and [3 = 10 -2 c), respectively. 
Other parameters are r = 10 3 , rj = 5 X 10 -2 and /i = 2.5 X 10 -2 , 
or equivalently, a = — 5 X 10 -2 , b± = 0 and a\ = 7.5 X 10 4 . The P2 
regime where k = 1 , induces the inversion of the initial condition 
at each round-trip as exemplified by the checkerboard pattern in 
d). It does not coarsen over longer time scales even for very large 
values of the asymmetry /3 = 0.1. The time trace is folded over 
times T in e) and also 2 T in f) for clarity. Parameters are 77 = 
—0.2 and n = —2.5 X 10 -2 , or equivalently a = —0.2, b\ = 0 and 
ai = 17.25 X 10 4 . The temporal profile is represented in g). Notice 
the large values of ai signaling that Eq. (3) remains qualitatively 
valid far from the bifurcation points. 


and defines, as in Eq. (2), whether or not the solution 
gets inverted from one round-trip to the next. For /3 = 0 
and a < 0, the solutions of the heteroclinic orbits of 
Eq. (3) read K± (. x) = ±\/H ptanh (|a| y/px ) where the 
± stands for the upward and downward domains walls. 
In deriving Eq. (1) we assumed that the symmetry break¬ 
ing term (3 scales as £ to enter as a perturbation of the 
solvability condition. Yet, although small, /3 has a deep 
impact upon the dynamics. 

We depict in Fig. la) an initial condition for Eq. (1) 
composed of an arbitrary succession of domains with val¬ 
ues corresponding to the plateaus of the PI or of the P2 
regimes. In the PI case (ft = — 1) as visible in Fig. lb,c), 
this multi-plateau pattern relaxes in a finite time to the 
upper (resp. lower) value when /3 < 0 (resp. /3 > 0 ). 
Such a coarsening scenario is very general and it is pre¬ 
served quite far from the bifurcation condition p = — 77 , 
and even when the discrete spectrum is also unstable, 
i.e. with p > 0. A multi-domain solution that verifies 
Eqs. (3,4) can be written as 

ip (x, n ) = K+ [x — Xq) + K_ (x — x± J (5) 

+ • • • + K_ (x - x~) - a/Kp 

with {x±} the ordered collection of the coordinates of 
the upward and downward kinks. For \/3\ <C 1, a varia¬ 
tional approach is justified [43], and we insert an ansatz 
as ip (x, s ) = K± [x — x± (s)]. By multiplying Eq. (3) by 


ip (x + 1 , n) = —Kip (x, n) 
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d x K± and integrating over the spatial coordinate, we find 
the effective equation for the motion of an isolated wall 
as 

^ —-— ( 6 ) 

dn 2 a^/\a\ 

demonstrating that the symmetry breaking term /? in¬ 
duces a splitting of the velocity of domain walls of op¬ 
posed “charge”, that will eventually collide. Notice that 
the value of the drift can also be found searching for 
heteroclinic solutions of Eq. (1) as in [35]. Both methods 
neglect the short range, exponentially decreasing interac¬ 
tions between nearby walls that is known to be attractive 
[16, 38, 43]. This can leads for (3 < 0 to a steady state, 
that is however unstable and cannot prevent the eventual 
collision and coarsening. Hence, one may conclude that 
it is impossible to store information in a scalar delayed 
system as any symmetry breaking nonlinearity leads to 
a coarsening of the information. 

We now demonstrate how a completely different regime 
can be obtained exploiting the anti-periodic solutions of 
Eq. (1) achieved by simply setting 77 < 0 to access the 
regime where k = 1 . We stress that here since fi < 0, 
bistability is lost and the only steady solution is T = 0 . 
Surprisingly, stable domains with well defined plateaus 
can be obtained even in this regime, although the time 
trace gets inverted every round-trip, thereby inducing a 
periodicity close to twice the delay value 2r, see Fig. ld- 
g). Although the solution inversion is interesting in its 
own right, such a behavior can be deduced intuitively 
from the singular mapping Eq. (2) setting k — 1. The 
striking result is that the domain walls remain motion¬ 
less, even for large (3 values, see Fig. Id) and Fig. le) 
where the checkerboard pattern visually disappears leav¬ 
ing only visible the transitions. We verified this robust¬ 
ness with other kind of symmetry breaking terms e.g. 
7 m T m with m = 0 and m = 4. One can still construct 
approximate analytic solutions of Eq. (3) where an up¬ 
ward or a downward kink must be complemented by an 
opposed kink at a spatial distance x = 1 in order to fulfill 
the anti-periodic boundary condition Eq. (4). A similar 
variational approach allows us to find the effective equa¬ 
tion of motion of an isolated wall, 

^ ~ (d x K+ + dxK-) (K+ + K-f dx = 0. (7) 


In other words, by inverting at each round-trip, the 
walls experience opposed drift velocities that cancel out 
explaining why similar results can be found for other kind 
of symmetry breaking nonlinearities. It was also demon¬ 
strated in [16] that the quadratic term cancels out of 
Eq. (3) even if it is large, see Eqs. (A. 17-18) in the ap¬ 
pendix A of [16] for details. Higher order nonlinearities 
in Eq. (1) potentially as large as ~ e 3_m m ay enters 
Eq. (3), but would cancel out similarly in Eq. (7). We 
conclude that P2 delayed systems are robust candidates 
for information storing as they are impervious to most 
experimental imperfections leading to a coarsening. 



Figure 2. (color online): a) Non linear output intensity / = \E\ 2 as 
a function of the injected field Y in blue. The black dots represents 
the amplitude of the upper and lower plateau in the P2 regime. 
The space-time diagram of the intensity in panel b) is folded with 
a period T = 2r + 0.1 and shows the stability of the domain walls, 
even in the presence of noise, see also panel c). Parameters are 
a = 2, J = —0.1, A = 2.3, fj = 0.1 and Q = 7 r. 

We study the case of an injected semiconductor laser 
below threshold subject to optical feedback, as domain 
stabilization in telecommunication lasers is very interest¬ 
ing also from an applicative point of view. Close, yet 
below, threshold and in the limit of weak optical feed¬ 
back, weak optical injection, small detuning and large 
delay, one can reduce the standard single-mode rate equa¬ 
tions of the semiconductor laser to the following delayed 
Ginzburg-Landau equation 

E = (1 + ia ) (j - |£| 2 ) E + iAE + Y + fie~ in E (t - t\ 8) 

see for instance [23] for details. Here, J denotes the de¬ 
viation from threshold of the bias current. The detuning 
between the frequency of the injection field ujy of ampli¬ 
tude Y (chosen real by definiteness) and the free running 
frequency of the laser scaled by the photon lifetime is 
A, whilst fj and £2 = c oyT are the amplitude and phase 
of the delayed optical feedback, respectively. At steady 
state (E = 0), the output power \E \ 2 as a function of Y 
can present a bistable S-shape response for some param¬ 
eters. We work close to the onset of bistability where the 
transition from low to high power has a sigmoid shape, 
see the blue line in Fig. 2a). 

Assuming that a P2 solution evolves between two 
plateaus whose values are E ±, it yields a system of equa¬ 
tions whose solution is represented by the black dots 
in Fig. 2a, b). Such a diagram suggests that these P2 
regimes are nascent from two saddle-node bifurcations 
of limit cycles and are also connected to the stationary 
solution (blue line) by two Andronov-Hopf bifurcations. 
However these bifurcations are subcritical meaning that 
the equivalent PDE would certainly take the form of a 
subcritical cubic-quintic Ginzburg-Landau equation. No¬ 
tice that for dimensional systems with a higher number 
of variables, as in Eq. (8), additional instabilities like the 
Eckhaus mechanism [44] may arise. Yet, Eq. (8) can be 
reduced to a single scalar delayed equation for the phase 
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Figure 3. (color online) (a) Experimental Setup. Legend: ML: 
Master Laser; PC: Polarization Controller; VOA: Variable Optical 
Attenuator; PM: Power Meter; FP: Fabry-Perot; CIRC: Circulator; 
PIN: 12 GHz Photodetector; SCOPE: 13 GHz Oscilloscope, (b) LI- 
curve of the SL at T = 293 K. (c) I/O power relationship of the SL 
under external injection with A/ = —10 GHz and / = 4.85 mA. 


in the limit of weak feedback and injection using averag¬ 
ing methods [45]. We show indeed in Fig. 2b, c) that this 
system is capable storing patterns of domain walls. 

The experimental setup is represented in Fig. 3. A 
1310 nm Fabry-Perot Laser (Slave Laser, SL) as used in 
telecommunication systems is subject simultaneously to 
external optical injection and delayed feedback after a 
round-trip time r = 65.4 ns. The optically injected sig¬ 
nal was generated with a tunable laser (Master Laser, 
ML). A polarization controller and a variable optical at¬ 
tenuator were respectively included after the ML to con¬ 
trol its polarization state and optical power level. Fig¬ 
ure 3b) plots the light-current (LI) curve of the solitary 
SL, showing a threshold current of I t h = 4.92 mA. Fig¬ 
ure 3c) depicts the SL’s input/output power relationship 
when subject solely to external injection. The device was 
biased below threshold with a current of I = 4.85 mA 
(i.e. I = 0.985 I t h) and an initial frequency detuning 
(A/ = finj—fFp) equal to —10 GHz was set between the 
frequencies of the injected signal and the resonance fre¬ 
quency of one of the SL modes. Figure 3c) illustrates the 
achievement of a gradual nonlinear switching response as 
the injection strength is increased from 0 to 31.5 /iW in 
agreement with the results of Fig. 2a). 

The SL was subsequently subjected to simultaneous 
optical injection from the ML (input power Pi n j = 
20 fi W) and delayed feedback. In this situation the sys¬ 
tem was operated in the sigmoid range as indicated by 
the blue square in Fig. 3c). The phase of the feedback 
can be changed by a small detuning of the injection field 
frequency. The large values of the time delay allows con¬ 
sidering these two parameters independent. 

Figure 4a) plots time trace over 5 round-trips (327ns). 



Figure 4. (color online). Temporal trace of the laser intensity a) 
and space-time diagram with folding parameter r b) and 2r c). 



Figure 5. (color online). Same as Fig. 4. A new pair of domain 
wall is nucleated at n = 128 in panel a) and remain stable while 
for different parameters bistable phase coarsening occurs. 


Figure 4a) shows that a temporal pattern which is in¬ 
verted every round-trip is obtained at the device’s out¬ 
put, in agreement with the predictions of Fig. lg) and 
Fig. 2c). The equivalent space-time diagrams for the 
time series of Fig. 4a) when the folding parameter is set 
approximately to r and 2 r are plotted respectively in 
Fig. 4b) and Fig. 4c) over a time window of 400 round- 
trips (i.e. ~ 26/is). The figures 4b-c) demonstrate 

the experimental achievement of anti-periodic dynamical 
regimes since peaks-troughs alternate every round-trip, 
see Fig. 4b). More importantly, the results also demon¬ 
strate the formation of stable domains of arbitrary size. 
The existence of various noise sources in the system in¬ 
duces a slow drift of the walls as seen in Fig. 4c), but 
not a coarsening if the walls remain at a sufficient large 
distance, in good agreement with theory. 

Noteworthy, we observed the nucleation of domain wall 
pairs after the arrival of an electrical perturbation into 
the system as shown in Fig. 5a) at n — 128. Additionally, 
a small change of the injection parameters allows to pass 
into the weakly bistable regime and to observe in Fig. 5b) 
the coarsening mechanism of the PI solutions, as in [35]. 
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In conclusion, we presented in this manuscript how a 
general property of delayed systems can be harnessed to 
prevent the domain coarsening in symmetry broken de¬ 
layed systems thereby allowing the storing of informa¬ 
tion. Such anti-periodicity has no equivalent in real spa¬ 
tially extended systems. We evidence the existence of 
stable domains in the coherent output of a semiconductor 
laser with optical feedback. These results offer exciting 
prospects for the controllable encoding of information. 
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